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Introduction 


•  Goal 

-  Develop  an  efficient  computational  strategy  to  predict  the  flame 
characteristics  and  thermal  response  of  aircraft  structures  immersed 
in  hydrocarbon  pool  fires 

-  Apply  strategy  to  crash  fire  scenarios  of  interest  to  the  FAA  and 
AFRL  (e,g.  past  accidental  fire  case  studies,  novel  commercial  / 
military  aircraft  platforms) 

-  Integrate  into  an  aircraft-crash-fire  simulation  framework 


■  Motivation 

-  Full-scale  aircraft  hydrocarbon  pool  fire 
testing  is  expensive  and  arduous 

-  The  FAA  is  concerned  with  unique  fire 
protection  challenges  New  Large  Aircraft 
(NLA)  pose  due  to  unusually  great 
dimensions,  fuel  quantities,  and  novel 
(composite)  materials 
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A  multi-tiered  simulation  framework  is  being  formulated  to  model  a  sequential  dynamic 
aircraft  crash  fire  event.  Applied  Research  Associates  SVO  is  developing  a  phase  I  aircraft 
crash  analysis  methodology,  with  the  AFRL/RXQD  following  suit  working  on  a  phase  II  liquid 
fuel  dispersion  and  phase  III  combustion  analysis  using  computational  fluid  dynamics. 
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Technical  Challenges 


»  Computational  efficiency 

-  Accuracy  vs.  Predictability 

*  Experimental  validation  data 

-  Confidence  Levels  /  Repeatability  /  Expense 

*  Multiple  non-linear  PDEs  governing  complex  flow  physics 

-  Combustion  /  Heat  Transfer  /  Turbulence  /  Multiphase  Flow 

*  Wide  range  of  spatial  and  temporal  scales 

-  Airbus  A380  9(1O0’s  m)  vs.  Turbulent  dissipative  length  scales  &(mm) 


HARDWARE  RESOURCES 


MULTI -SCALE 
TURBULENT  FIRE  PLUME 


MESH  TOPOLOGY 
OPTIMIZATION 
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Progress  Path 


•  CFD  Model  Development 

-  Computational  Resources  /  Physics  /  Boundary  Conditions  i  Domain 
Development 

•  Model  Validation 

-  J.M.  Suo-Antilla  and  L.A.  Gritzo.  “Thermal  Measurements 
from  a  Series  of  Tests  with  Large  Cylinder  Calorimeter  on  the 
Leeward  Edge  of  a  JP-8  Pool  Fire  in  Cross-Flow.”  SAND 
2001-1986. 

-  Low  /  Medium  /  High  cross  wind  cases 

•  Application  to  Full-Scale  Aircraft 

-  Boeing  707 

-  NLA  {Airbus  A380,  Boeing  707,  Tyndall  AFB  Mock-Up,  etc.) 

-  Future  Aicraft 

•  Integration  into  a  Dynamic  Crash  Aircraft-Crash-Fire 
^  ’  Analysis 


The  bullets  are  faded  at  the  bottom  of  slide  5  because  it's  a  progress  path  and  those 
items  have  not  been  completed  yet. 
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ANSYS  Fluent  12.1:  11  Equation  Model 

-  3-D  Navier-Sfokes  (5) 

-  Non-Premixed  Combustion  (Mixture  Fraction  PDF)  Approach  (2) 

-  Realizable  k-e  (RANS)  Turbulence  (2) 

-  Discrete  Ordinates  (DO)  (1 ) 

-  Single  Step  Khan  and  Greeves  Soot  Model  (1 ) 


PROBABILITY 
DENSITY  FUNCTION 
(PDF) 

ANSYS  Fluent 
Theory  Guide 
(2008) 


TURBULENT  - 
CHEMISTRY 
INTERACTION 


DO  RADIATION 
DISCRETIZATION 


llCFD  in  Fire 
Engineering/1  Ye  oh 
and  Yuen  (2008) 


CFD  Model  Validation:  Experiment 


•  Low  /  Med  /  High  Cross  Winds  Over  a  Cylinder 


D  AQ  LAYOUT  E  XPER I  M  ENTAL  SET-UP  MEDIUM  CROSS-WIND  (INSTANT  AN  ECUS) 


Case  Name 

Time  After  Ignition  (s) 

Wind  Velocity  (m  s'") 

Wind  Direction  O 

Low 

225  -  350 

1.9  ±0.2 

-36.9  ±  5.7* 

Medium 

400-575 

5.4  ±  1.2 

-1 1.4  ±  12.5* 

High 

300-600 

10,2  ±1,7 

-22.7  ±  B.3® 

Sandia  National  labs  has  developed  an  experimental  test  series  designed  specifically  to 
provide  flame  and  skin  temperature  and  heat  flux  data  on  a  cylinder  immersed  in  a  jet  fuel 
fire  under  varied  wind  conditions.  Certain  segments  of  their  data  was  extracted  and 
averaged  due  to  optimal  quasi-steady  wind  conditions.  A  low,  medium,  and  high  wind  data 
set  was  selected  from  their  data  reduction  to  formulate  a  comparable  modeling 
environment  to  compare  results. 
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CFD  Model  Validation:  Model  Set-Up 


|  Clt£9Sw»E  OUTLET 

|  INLET 

VtOUTOFPAPEftg 

NEAR  FIELD 
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CROSS  WIND 
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Fire  Inlet  Condition 

Value 

Flow  Rate  (kgitrV) 

0.058 

Temperature  (K) 

510 

1±±jLUl±jLL 


JET  A  VAPOR  IN  FLOW 
(MASS  F  LOW  INLET  Q  510  K) 


VZ/77777Z77777Z7777/ 

6ROURO 
(WALL®3MK) 


Case  Name 

Mean  Wind  Velocity  (ms1) 

TJ  (%) 

Reynolds  Number  (Ke^) 

Low 

1,9 

10.5 

4.76  x  10* 

Medium 

5.4 

22.2 

1.35  x  10* 

High 

10.2 

16.7 

2.56  x  10* 

N'i>ti :  Thti  £hafatt*rialii  (urhultnl  Irnyth  itak  uatid  far  All  caar»  is  0.2M  m 


The  following  boundary  conditions  were  set  reflecting  measurements  and  observations 
made  by  Sandia  National  labs  to  create  similar  physical  conditions.  Far  field  model 
boundaries  were  defined  to  simplify  model  orientation  with  respect  to  the  oncoming 
measured  wind  direction  from  experiments.  The  fire  inlet  conditions  were  extrapolated 
from  Sandia  experimental  data  and  turbulent  boundaries  conditioned  to  reflect  the 
fluctuations  measured. 

All  wall  conditions  were  initialized  to  300  K,  but  were  allowed  to  float  (or  rise)  during  the 
CFD  solution  process  to  account  for  combustion  convective  and  radiant  heat  transfer 
effects. 

Reynolds  number  calculations  are  based  upon  the  cylinder  diameter,  and  reflect  an 
increasing  turbulent  flow  regime  due  to  increasing  mean  velocity  validating  the  need  for 
turbulent  versus  laminar  flow  modeling  assumptions. 
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Pointwise  Gridgen  15:  Multi-Block  Hybrid  Topology 

-  Structured  (hexahedral)  high  aspect  ratio  cells  used  for  far  field 
atmospheric  boundary  layer  development  and  cylinder  boundary  layer 
growth 

-  Unstructured  (tetrahedral)  cells  used  to  link  structured  blocks 


-  400K  NODES 


CFD  Model  Validation:  Model  Set-Up 
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Fire  Intel  Condition 

Value 

Flow  Rale  (kg' nr V) 

0.058 

Temperature  (K) 

510 

JET  A  VAPOR  IN  FLOW 
(MASS  F  LOW  INLET  Q  510  K) 


77/77777//77777//77777s 

GROUND 
[WALL  ®  304  K) 


Case  Name 

Mean  Wind  Velocity  (ms1) 

n  m 

Reynolds  Number  (Ke^i) 

Low 

1,9 

10,5 

4.76  x  10* 

Medium 

5,4 

22,2 

135  x  10* 

High 

10.2 

16,7 

2.56  x  10* 

Note:  Th«  £hafatt*rialii  (urhultnl  Irnyth  dcata  UTiud  lor  all  caar»  is  0.2M  m 


Low  wind  flame  temperature  measurements  compare  very  well  with  model  results, 
especially  when  taking  into  account  dissimilarities  between  experimental  set-up  and  model 
conditions. 

The  low  cross-wind  conditions  create  the  best  conditions  for  a  pure  diffusion  flame 
compared  to  the  higher  wind  conditions,  reflecting  an  environment  most  representative  of 
the  model  design  assumptions. 
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The  medium  wind  comparative  condition  compares  reasonably  well  like  the  low  wind 
condition,  with  the  most  normally  distributed  wind  gust  across  the  cylinder.  With  the 
highest  turbulence  intensity  measured  of  all  three  cases,  this  case  provides  the  most  even 
and  well-averaged  results  most  favorable  for  turbulent  RANS  modeling  conditions. 
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f  CFD  Model  Validation:  High  Wind 


Flame  Temperatures 


ItJ  MS'Mi.r 

iJt'Kl 


io,iw&V-ai.r 
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The  high  wind  validation  case  shows  the  highest  discrepancy  of  all  three  cases  for  2  primary 
reasons.  Here,  the  convective  force  of  the  cross  wind  condition  begins  to  compete  heavily 
with  the  diffusion  combustion  process  creating  a  fire  plume  nearing  the  modeling 
assumption  limitations.  In  addition,  increased  cross  winds  create  greater  flow  separation 
around  the  leeward  side  of  the  cylinder.  With  the  employment  of  a  reduced  order  (2- 
equation  linear  eddy-viscosity)  turbulence  model  to  save  computational  cost,  cylinder  flow 
separation  and  ultimate  detachment  becomes  increasingly  hard  to  predict,  shifting  average 
flamelet  location  off  of  the  cylinder  surface.  This  creates  large  discrepancies  between  skin 
temperatures  and  heat  fluxes  leading  to  the  largest  divergence  between  experimental  and 
model  results. 
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CFD  Model  Validation:  Low  Wind 


•  Skin  Temperatures 


SUO-ANTILLA  &  GRITZO  (2001)  -  EXP  CFD 


Low  speed  skin  temperature  modeling  results  reflect  reasonably  well  compared  to 
experiment  once  again  when  taking  into  account  all  of  the  model  simplifications.  Maximum 
temperature  magnitudes  are  achieved  and  trends  are  generally  comparable. 
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Similar  to  the  low  wind  case,  medium  wind  conditions  create  favorable  results  as  well  with 
applied  conditions  still  well  within  model  assumptions. 
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CFD  Model  Validation:  High  Wind 


Skin  Temperatures 
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Model  issues  with  flow  separation  become  largely  apparent  here  with  windward  surface 
temperatures  comparing  reasonably  well  with  experimental  results,  but  the  leeward 
cylinder  side  diverging  significantly  due  to  the  flame  lifting  off  of  the  surface. 

Flame  lifting  can  be  caused  by  a  combination  of  both  turbulent  flow  model  breakdown 
combined  with  combustion  model  assumption  divergence. 
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Heat  flux  magnitudes  tend  to  follow  suit  with  skin  temperature  measurements,  with  the 
low  wind  condition  creating  acceptable  heat  flux  magnitudes  and  trending  profiles 
between  model  and  experiment. 
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Although  heat  flux  model  and  experimental  magnitudes  aren't  as  comparable  compared  to 
the  low  wind  speed  case,  the  profiles  compare  well  showing  major  flow  structures  are  well 
captured. 
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Leeward  heat  flux  comparisons  between  model  and  experiment  suffer  the  most  for  the 
high  wind  regime  for  the  same  reasons  discussed  previously.  Heat  flux  magnitudes  and 
trending  profiles,  however,  due  compare  well  on  the  windward  side  where  flow  separation 
is  not  an  issue. 
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The  following  three  pictures  depict  the  average  flame  shape  and  temperature  in  the  3 
cross-cuts  compared  against  experimental  data,  in  addition  to  cylinder  skin  temperatures. 

The  low  wind  case  shows  the  lowest  plume  temperatures  largely  due  to  the  least  amount 
of  fuel  air  mixing  combined  with  the  largest  off  normal  wind  direction  with  respect  to  the 
cylinder  surface. 

The  medium  wind  case  has  the  highest  fluctuating  wind  condition  and  the  most  normal 
wind  direction  creating  the  most  definitive  plume  shape  and  increased  temperatures. 

The  high  wind  case  had  the  least  amount  of  wind  fluctuations  but  the  highest  magnitudes 
causing  increased  swirl  behind  the  cylinder  and  non-coherent  plume  structures. 
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CFD  Model  Application:  Aircraft 


•  Full-scale  Boeing  707  selected  as  baseline  aircraft 

•  NFPA  403  PCA  fire  condition  under  low  speed  wind  conditions 


AIRCRAFT  SHAPE/SIZE  COMPARISON  NFPA  403  FIRE  CONDITION 


The  Boeing  707  was  selected  as  the  baseline  aircraft  due  to  its  comparable  size  and  shape 
with  the  average  commercial  jetliner. 

The  PCA  fire  condition  as  described  in  NFPA  403  was  used  for  model  conditions,  with 
formal  boundary  conditions  derived  from  the  low  to  medium  cross-wind  speeds  that 
compared  the  best  in  the  validation  study. 
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Flow  visualization  of  the  Boeing  707  in  a  5  mph  (equivalent  low  speed  condition)  with  10% 
turbulence  intensity  is  depicted. 

Similar  to  the  validation  cylinder  results,  thermal  surface  conditions  are  affected  by  the 
local  flow  structures  created  by  the  turbulent  mixing  between  the  combusting  fuel  vapors 
and  incoming  oxygenated  wind.  More  complex  surface  shapes  create  a  more  varied  profile 
creating  localized  hot  and  cold  spots  due  to  varied  aerodynamic  shapes  (stream-wise  main 
wings  versus  the  bluff  body  the  vertical  tail  presents  to  the  oncoming  cross-flow). 
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Conclusions 


*  Experimental  flame  and  skin  temperatures  compare  well  profile 
trend  and  magnitude-wise  to  CFD  at  low  to  medium  wind 
velocities,  with  the  largest  discrepancies  in  the  high  wind  regime 

*  Cylinder  flow  separation  appears  to  be  a  first  order  turbulent  flow 
modeling  effect  leading  to  the  divergence  on  the  leeward  cylinder 
surface  thermal  effects 

*  Sharp  flamelet  gradients  are  indicative  of  fast  chemical  kinetic 
time  scales  largely  absent  in  the  “mixed-is-burnt”  model  based 
upon  slower  fluid  motion  time  scales 

*  Simple  soot  modeling  provides  sufficient  radiant  heat  absorption 
to  accurately  meter  total  heat  transfer  magnitudes,  but  lacks  in 
qualitative  plume  size  and  shape  prediction 

*  Full-scale  aircraft  thermal  and  fluid  structures  trend  similarly  to 
the  validation  study,  with  resource  needs  increasing  by  about  4x. 
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Future  Work 


•  Analyze  alternative  CFD  sub-models  (combustion  /  turbulence  / 
soot)  to  increase  accuracy 

•  Gather  more  validation  I  material  property  data 

•  Run  unsteady  cases  to  determine  transient  effects  and  examine 
material  heat  up  times 

-  Complete  Airbus  A380  and  Boeing  787  case  studies  (Possibly 
Tyndall  AFB  Airbus  Mock-Up?) 

•  Apply  to  alternate  aircraft  I  combustion  scenarios  of  interest 

•  Ultimate  integration  into  an  aircraft-crash-fire  simulation 
framework. 
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